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Abstract 

An algorithmic method using conservation law multipliers is introduced that yields 
necessary and sufficient conditions to find invertible mappings of a given nonlinear 
PDE to some linear PDE and to construct such a mapping when it exists. Previous 
methods yielded such conditions from admitted point or contact symmetries of the 
nonlinear PDE. Through examples, these two linearization approaches are contrasted. 
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1 Introduction 

In previous work [1-3], algorithmic methods were presented to determine whether or not 
there exists an invertible mapping of a given nonlinear partial differential equation (PDE) 
to some linear PDE through computing the point or contact symmetries admitted by the 
nonlinear PDE. Moreover, it was also shown how to construct such an invertible mapping, 
when one exists, from the admitted symmetries. 

In this paper, we present an alternative algorithmic method for mapping a nonlinear 
system of PDEs invertibly to some linear system of PDEs by using admitted conservation 
laws of the nonlinear system. Any linear system possesses an infinite set of conservation law 
multipliers satisfying its adjoint system. This feature of linear systems will be exploited to 
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detect whether or not a given nonlinear system can be linearized by an invertible transfor- 
mation and, when such a linearization mapping exists, to obtain the explicit form of the 
linearizing transformation from an associated conservation law identity. 

In particular, using the determining equations for conservation law multipliers, one can 
find whether or not an invertible linearization mapping will exist, and also find the adjoint 
of a target linear system as well as the transformation of the independent variables in the 
mapping whenever an invertible linearization is possible [4]. Most importantly, the new work 
presented in this paper will show that a conservation law identity coming from the multiplier 
equations of an augmented system consisting of the given nonlinear system and the adjoint 
of the linear target system can be used to determine the transformation of the dependent 
variables (in addition to that of the independent variables) in the invertible linearization 
mapping explicitly. This work thus allows one to detect and construct linearization mappings 
completely by the use of algorithmic methods for obtaining multipliers of conservation laws 
[5-9]. 

Our results generalize earlier computational work [10] in which a fully automatic algo- 
rithm was given using the computer algebra program Crack/Conlaw [11] to detect and 
construct linearizations in the case when the mapping only transforms the dependent vari- 
ables of the given nonlinear PDE system. The work in [10] also gave an extension of this 
algorithm to find partial linearizations that consist of factoring out a linear differential op- 
erator from all PDEs in a given system (with the operator being of lower order than the 
highest order of derivatives in the system). 

To begin we summarize two well-known theorems characterizing invertible mappings. 
Our notation follows that in [12], with the conventions that Latin indices i,j, k will run 1 
to n (labeling independent variables), and Greek indices u, /i, A and a, (3, 7, a, r will run 1 to 
M and 1 to m (labeling PDEs and dependent variables), respectively. Throughout we use 
upper/lower case notation to distinguish arbitrary functions U, V from functions u,v given 
by solutions of a PDE system. 

Let H{x, u} denote a given /cth-order nonlinear system of M PDEs with n independent 
variables ) and m dependent variables u = (u 1 , . . . , u m ). In particular, let the 

PDEs in R{x, u] be given by 

G"[u] = G u (x,u,d x u,...,d^u) = 0, v = 1, . . . , M, (1.1) 

where d x u denotes jth-order partial derivatives of u with respect to x. 

Let S{z,w} denote a /cth-order linear target system of M PDEs with n independent 
variables z — (zi, . . . , z n ) and m dependent variables w = (w 1 , . . . , w m ), 

H u [w} = H u (z,w,d z w,...,d*w) = 0, u=l,...,M, (1.2) 

where d{w denotes jth-order partial derivatives of w with respect to z. 

Note that for common PDE systems such as heat or diffusion equations and wave equa- 
tions, as well as more general parabolic, hyperbolic and Hamiltonian equations, and elliptic 
equations, we will have m = M. However, for complete generality we allow m 7^ M hereafter, 
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as would occur for instance if one considers a PDE system of dynamical evolution equations 
with differential constraints. 

Now if one seeks an invertible mapping \i of a given nonlinear PDE system f 1 1.1 1) to a linear 
target PDE system (jl.2p . then the following two theorems show that the mapping is restricted 
to being a contact transformation when (11. ip and (jl.2p both contain a single dependent 
variable (m = 1; which we will call the scalar case) or to a point transformation when fll.il) 
and (11.21) each contain two or more dependent variables (m > 2; the multicomponent case). 

Theorem 1.1 (Backhand [13]). If u = u 1 is a scalar (m = 1), then a mapping \i defines 
an invertible transformation from (x, u, d x u, . . . , d^u) - space to (z, w, d z w, . . . , d^w) - space 
for any fixed p if and only if fi is a one-to-one contact transformation of the form z = 
4>(x,u,d x u), w = ip(x,u,d x ii), d z w = p(x,u,d x u), which is subject to the contact condition 
Dxip = pD x <P where D x denotes total derivatives with respect to x. 

Theorem 1.2 (Miiller and Matschat [14]). If u = (it 1 , . . . ,u m ) has two or more components 
(m > 2), then a mapping fi defines an invertible transformation from (x,u,d x u, . . . ,d%u) - 
space to (z, w, d z w, . . . , d p z w) - space for any fixed p if and only if jj, is a one-to-one point 
transformation of the form z = <f)(x, u), w = ip{x, u). 

For formulating necessary and sufficient conditions for the existence of such mappings 
between PDE systems, we will consider only the natural situation where the given nonlinear 
system (II. ip and the linear target system (ll.2p contain the same number M of PDEs with 
the same differential order k, and no further assumptions will be required (such as local 
solvability, or involutivity). But for the computational implementation of our results it is 
important that the given system of PDEs (II. ip be consistent in the sense that all integra- 
bility conditions are satisfied modulo these PDEs. (To make this precise, one would have 
to prescribe an elimination algorithm which requires the definition of a total ordering of 
derivatives and of monomials built from them, as discussed in [10].) 

The rest of this paper is organized as follows. In section two, we review the method 
presented in [1-3] to linearize nonlinear PDEs through admitted symmetries. In section 
three, we present the new method to linearize nonlinear PDEs through admitted conservation 
law multipliers and outline the computational steps involved in this method. As examples to 
illustrate and contrast both methods, in section four we consider linearizations of Burgers' 
equation, a pipeline flow equation, and a nonlinear telegraph system. We make some closing 
remarks in section five. 

2 Use of Symmetries to Construct Linearizations 

Necessary and sufficient symmetry conditions are stated in [1-3] to determine whether or 
not a given nonlinear system R{x, u} of PDEs (11.11) can be transformed to some linear 
target system S{z, w} of PDEs (ll.2p by an invertible mapping \i. These conditions provide 
a symmetry-based algorithmic procedure to determine if they hold for a given nonlinear 
system of PDEs and (when these conditions do hold) a second symmetry-based algorithmic 
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procedure to construct such an invertible mapping /i. To apply these methods to a given 
nonlinear system of PDEs, it is unnecessary to know a specific linear target system of PDEs. 
In particular, the linear target system (when one exists) will arise from computing the 
admitted point or contact symmetries (Lie group of point or contact transformations) of 
the given nonlinear system (first procedure). Moreover, these admitted symmetries yield 
equations to construct a specific mapping /i (second procedure). 

The starting point is that a linear system S{z, w} of PDEs, given by 

L[z}w = g(z) (2.1) 

as defined here in terms of a linear operator L[z], with any inhomogeneous term g(z), is 
completely characterized by its admitted infinite-parameter set of point symmetries 

X = /W^, (2-2) 
where f(z) is any function satisfying the related linear homogeneous system 

L[z)f = 0. (2.3) 

A point transformation is mapped to another point transformation under any specific 
point transformation (x,u) — > (z,w), and a contact transformation is mapped to another 
contact transformation (which could be a point transformation) under any specific contact 
transformation (x,u,d x u) — > (z,w,d z w) (which could be a point transformation). Hence an 
invertible point transformation mapping a nonlinear system R{x, u} of PDEs to some linear 
system S{z, w} must map each point symmetry admitted by R{x, u} to a point symmetry 
admitted by S{z,w}; an invertible contact transformation mapping a nonlinear scalar PDE 
system R{x, u] to some linear scalar PDE system S{z, w}, must map each contact symmetry 
admitted by R{x, u} to a contact symmetry (which could be a point symmetry) admitted 
by S{z,w}. Moreover, this mapping must be an isomorphism, i.e., there is a one-to-one 
correspondence between the point (contact) symmetries admitted by H{x,u} and those 
admitted by S{z,w}. Consequently, if there exists an invertible transformation that maps 
the nonlinear PDE system R{x,m} to some linear PDE system S{z,w}, then R{x,m} must 
necessarily admit an infinite set of point (contact) symmetries. 

The following two theorems on mappings of nonlinear PDEs to linear PDEs were pre- 
sented and proved in [1-3]. 

Theorem 2.1 (Necessary conditions for the existence of an invertible linearization mapping). 

If there exists an invertible mapping fx of a given nonlinear system R{x, u} of PDEs (11. ip to 
some linear system S{z,w} of PDEs (11.21) . then in the multicomponent case (i.e. m > 2): 

1. Li is a point transformation of the form 

Zi = <f>i(x,u), i = l,...,n, w a = 4> a (x, u), <r = l,...,m; (2.4) 



4 



2. R{x, u] admits an infinite-parameter Lie group of point symmetries given by an in- 
finitesimal generator 



X = Ux ' u) lk +r]T{x > u) ^ (2 ' 5) 
with infinitesimals of the form 

& = a i(T {x,u)F°{x,u), rf = p T a (x,u)F' 7 (x,u), (2.6) 

where ai a , (3^ (i = 1, . . . , n; a, r = 1, . . . , m) are specific functions of x and u, and 
where F = (F 1 , . . . , F m ) is an arbitrary solution of some linear system of PDEs 

L[X]F = (2.7) 

in terms of some linear differential operator L[X] and specific independent variables 
X = (Xi(ar,it), • • .,X n (x,u)) = (fa, . . .,(p n ). 

Theorem 2.2 (Sufficient conditions for the existence of an invertible linearization map- 
ping). Suppose a given nonlinear system H{x,u} of PDEs U.l\) admits infinitesimal point 
symmetries (I2.5P of the form (12.61) involving an arbitrary solution F(X) of a linear system 
( 12. 71) with specific independent variables X = (X\{x,u), . . . ,X n (x, u)). If the m first order 
scalar linear homogeneous PDEs 

Oi <r (x,u)^- + /^(a;,u)^: = 0, a = 1, . . . ,m, (2.8) 

whose coefficients are formed from (12.61) have <f>\ = X\(x, u), . . . , <fr n = X n (x, u) as n func- 
tionally independent solutions, and if the m 2 first order linear inhomogeneous PDEs 

dip T dip T 

aia ( x,u ^~dx~ + ^ X ' u ^lhn = ^' T,a = lj ' ' ' ' m ' ^ 2 ' 9 ^ 

(where 5^ is the Kronecker symbol) have a particular solution ip = (ip l (x,u), . . . ,ip m (x,u)), 
then the mapping /i defined by 

Zi = Xi(x,u), i = l,...,n, w a = ip a (x,u), a = l,...,m, (2-10) 

is invertible and transforms H{x,u} to the linear system S{z, w} of PDEs given by L[z]w = 
g(z), for some inhomogeneous term g(z). 

It is easy to see that Theorems 2.1 and 2.2 also hold in the scalar case (i.e. when 
m = 1) if one seeks to transform nonlinear scalar PDEs to linear scalar PDEs only by a 
point transformation. Furthermore, these theorems in the scalar case can be extended to the 
more general possibility of seeking invertible mappings of nonlinear scalar PDEs to linear 
scalar PDEs by contact transformations. For details, see [1,2]. 
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3 Use of Conservation Law Multipliers to Construct 
Linearizations 



Now we present a new method for finding mappings that transform a nonlinear system of 
PDEs invertibly to some linear system of PDEs. Our method uses admitted conservation law 
multipliers of the nonlinear system and is based on two theorems that are the conservation 
law multiplier-based counterparts of the symmetry-based Theorems 2.1 and 2.2. 

The first new theorem gives necessary conditions that conservation law multipliers must 
satisfy in order for a given nonlinear system R{x, u] of PDEs (11.11) to have an invertible 
mapping to some linear target system S{z, w} of PDEs (11. 2p by a point transformation. As 
in the case of Theorem 2.1, the necessary conditions yield an algorithmic procedure to find 
a specific linear target system as well as its independent variables. 

The second new theorem gives sufficient conditions for the existence of the mapping. In 
particular, as in the case of Theorem 2.2, these conditions yield an algorithmic procedure to 
find the components of an explicit point transformation giving an invertible mapping to the 
linear target system. 

The extension of these two theorems to include contact transformations will be discussed 
in section 3.2. 

3.1 Necessary and Sufficient Linearization Conditions 

Recall, a conservation law of a system R{x, u] of PDEs (11.11) is a divergence expression 
DTi[u]/Dxi that is required to vanish for all solutions u(x) = (u (x), . . . , u m (x)) of the 
given PDE system. Conservation laws are nontrivial if and only if the conserved fluxes Tj [u] 
do not have the form of curls DQij[u}/ Dxj, with Oij[w] = —Qji[u], for every solution u(x). 

Definition 3.1. Conservation law multipliers for a system H{x,u} of PDEs (II. ip are a 
set of factors {A u (x, U, d x ll, . . . , d l x U)} [y = 1,...,M) such that, for arbitrary functions 



holds as an identity for some functions T = (Tift/], . . . , Y„ [£/]), where each factor A„[L/] is 
non-singular for all solutions U = u(x) of the given PDE system (11.11) . 

It is well known that multipliers will yield all nontrivial conservation laws (to within the 
addition of curls to the conserved fluxes) admitted by a given PDE system R{x, u] if it is 
of Cauchy-Kovalevskaya type. 

For any PDE system H{x,u}, its admitted conservation law multipliers satisfy the de- 
termining equations 



U= (U 1 (x),...,U m (x)) 



k v [U]G v \U] = DTi[U]/Dxi 



(3.1) 



E u „(A v [U]G v [U]) = 0, a = l,...,m, 



(3.2) 
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given by the Euler operators 



d D d D 2 d , _ t D k d 

'->u° 



En« = ^7Z- 77T7 + ^ ^ " +■■■ + (~l) k ^ ^ 7777~~ + 



dU° DxidUZ. DxiDxjdU^ Xj "' v ' Dx { . . . D Xj dU°. 

which annihilate divergence expressions. One can solve these determining equations by an 
algorithmic procedure [8, 10] analogous to that for solving the determining equations for 
point or contact symmetries in evolutionary form. (The similarity is made explicit in the 
work in [6, 7] which shows how to reformulate (13.21) as a kind of adjoint of the symmetry 
determining equations for Cauchy-Kovalevskaya PDE systems.) 

The starting point, leading to the formulation of an invertible mapping between nonlinear 
and linear PDE systems through conservation law multipliers, is that any linear operator 
L and its adjoint operator L* satisfy a conservation law identity. Specifically, consider 
a fcth-order linear PDE system Lfzjty = 0, denoted by S{z, w} with dependent variables 
w = (w l , . . . ,w m ) and independent variables z = (z l , . . . , z n ). In particular, let the linear 
PDEs be given by 

L u a [z]w a = 0, z/ = l,...,M (3.3) 
in terms of linear operators 

d d k 

L v a [z\ = b va (z) + b vai (z) — H h b vail ...i Az)- — . (3.4) 

dzi dz h ■ ■ ■ dz ik 

The corresponding adjoint linear system, L*[z]t> = 0, is given by 

LT[zV = W*K ~ ^(W*K) +■■■ + (-!) fc ^. ...g (W--^K) (3 5) 
= 0, a = 1 m, 



where L* a [z] is the linear operator defined through formal integration by parts of the oper- 
ator L u a [z] with respect to z. Then for arbitrary functions V(z) = (V 1 (z), . . . , V M (z)) and 
W(z) = (W 1 (z) 1 . . . , W m (z)) ) which can be viewed as multipliers {V^(z)} and {VT Q (z)} for 
the augmented linear system consisting of the linear system (13.31) and the adjoint system 
(13. 5p . there is a conservation law identity 

DTi/Dzi = 5 u ,V»U a [z\W a - 8 <xP W p Ll a [z]V v (3.6) 



holding for some specific functions YjfV 7 ,!^] that have a bilinear dependence on the multi- 
pliers and their derivatives. (An explicit expression for the fluxes Tj[V, W] is given in [15].) 
Here both 5 UfJi and 5 a p are Kronecker symbols. 

a )We freely raise and lower indices for convenience here in displaying the adjoint relation between L and 
L*. 
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Propostion 3.1. Suppose a given nonlinear system R{x, u] of PDEs Al.l\) can be mapped 
into some linear system S{z,w} of PDEs A3. 3\) by some point transformation \2. 4\ ). Then 
the nonlinear and linear PDEs will be related through a set of factors {Q^(x,U,d x U)} such 
that 

Qt[U]G"[U)=I£[z]W a , fi=l,...,M, (3.7) 

for arbitrary functions U(x) = (U l (x), . . . , U m (x)), where the functions W(z) = 
{W 1 {z) ) . . . , W m (z)) are given by the linearization mapping 

z = (j)(x,U(x)), W(z) =ip(x,U(x)). 

(Explicit expressions for can be readily derived from the mapping formulae given in [16].) 
Moreover : 

1. This mapping will define an invertible point transformation iff {<f>i(x,U(x))} (i = 
1, . . . ,n) are functionally independent and {Q„(x, U(x), d x U(x)) (v,ix— 1, . . . , M) are 
non- degenerate, namely, det(D<f)i(x,U)/Dxj) ^ and also det(Q£(x,U,d x U)) ^ 0. 

2. Under such a point transformation the conservation law identity A3. 6\) becomes 

DVi/Dxi = (5 X ,V X Q>[:[U]G U [U] ~ 6 a/3 W p L* u a [z}V u )J[U], (3.8) 
where DTi/Dxi is related to DTi/Dzi by the non-vanishing Jacobian factor 

(x,U(x)) 



J(x,U,d x U) 



Dz 




= det ^ 


Dx 





Dxj 



(3.9) 



(See [12] for the explicit expression for Tj in terms o/TjJ 



This proposition leads to the following two theorems for linearization from admitted con- 
servation law multipliers, which are the counterparts to Theorems 2.1 and 2.2 for linearization 
through admitted point symmetries. 

Theorem 3.1 (Necessary conditions for the existence of an invertible linearization mapping). 

If there exists an invertible point transformation (12. 4p for m > 1 under which a given kth- 
order nonlinear system R{x, u] of PDEs (11.11) is mapped to some linear system S{z, w} of 
PDEs A3. 3\) . then R{x, u} must admit conservation law multipliers of the form 

K[U] = 5 x y(X)QZ[U] J[U], (3.10) 

where f/i, v = 1, . . . , M) are specific functions of x, U and d x U; v = (v 1 , . . . , v M ) is an 
arbitrary solution of some linear system of PDEs 

L;[I]/ = ^(I)/ + & a(i! (I)4 i + . . . + 6 a/lii ... it (I)^ ...^ = 0, a = l,...,m, (3.11) 

in terms of specific independent variables X = (Xi(x, U), . . . , X n (x, U)); and J[U] = 
\DX/Dx\ = det(DXi(x,U)/Dxj) which is a function of x, U and d x U . 
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Proof. The existence of multipliers (13.101) follows from equation (13.81) with X playing the 
role of z, L^[X] playing the role of L* u a [z\ and with = satisfying system (13. lip . □ 

Theorem 3.2 (Sufficient conditions for the existence of an invertible linearization mapping). 
Suppose a given nonlinear system R{a;, u} of PDEs (11.11) admits conservation law multipliers 
of the form (13.101) where the components of v are dependent variables of some linear system 
(13. lip with specific independent variables X = (Xi(x,u), . . . ,X n (x,u)). Let L*[X) be the 
adjoint of the linear operator L[X] in Ii3.11\) and consider the augmented system of PDEs 
consisting of the given nonlinear system (11.11) and the linear system ( 13.111) . Then, from 
identity ( 13. 8p . there exists an infinite set of multipliers 

{A„[U,V],A a [U\} = {Ql(x,U,d x U)5 Xv V x (X(x, U))J[U], -6 af} W p (x, U)J[U}} 
(/i = 1, . . . , M; a — 1, ... ,m) yielding a conservation law 

A V [U, V]G V [U] - A a [U]L«[X(x, U)]V» = DT l /Dx l (3.12) 
for some specific functions Ti, in terms of the Jacobian 

J[U] = \DX/Dx\ = det(DXi(x, U)/Dxj). (3.13) 
This conservation law ( I3.12p is equivalent to the identity 

5 x ,V x QZ(x, U, d x U)G»[U] - 5 aP W a (x, U)L?[X(x, U)]W = DT./DX,, (3.14) 

holding for some functions Tj. // the variables {X{\ (i = l,...,n) are functionally in- 
dependent and if the factors {Q„} (\i,v = 1,...,M) are non-degenerate, then the point 
transformation given by 

z = X(x, u), w = W{x, u) (3.15) 

maps the nonlinear system of PDEs (II .ip invertibly into the linear system S{z,w} given by 

L*:[z]w<* = b au {z)w a - —(b avi {z)™°) +■■■ + (~l) fc ^ ...^ (b auil ... ik (z)w a ) (3 lg) 
= 0, u = l,...,M. 

Proof. Since L[X] is a linear operator, the identity (13. 6p gives 

S aP W a (x, U)i%[X(x, U)]V" = 5^V»L* a "[X(x, U)}W a (x, U) + DVi/DXi (3.17) 
for some specific functions Tj [U, V, W] . Consequently, the identity (13.141) becomes 

5^V\Qt\U}G v \U] - LT[X(x, U)]W a (x, U)) = D(T t + T t )/DX t . (3.18) 
Now apply the Euler operators with respect to i.e., 

d D d 
Evv = t; t; — tt + • • • i /i = l,..., M, 

v dW DXi dVg. ' P ' ' ' 



to each side of equation (I3.18P . Each of these Euler operators annihilates the right-hand side 
of (13.181) and hence we obtain the identity 



Qt[U]G v [U\ = L?[X{x, U)]W a (x } U) (3.19) 

holding for arbitrary functions U. Now suppose U = u(x) solves the given nonlinear system 
of PDEs (II. II) . then it follows that w = W(x, u(x)) solves the linear system given by (13. 16ft . 
Consequently, one obtains the point transformation (13.151) which maps the nonlinear system 
of PDEs (II. ip into the linear system (13.161) . Invertibility of this transformation is established 
by Proposition 3.1. □ 



3.2 Computational Steps 

We now outline the computational steps involved in applying Theorems 3.1 and 3.2 to 
linearize a given nonlinear system of PDEs through admitted conservation law multipliers. 

Step 1: For a given nonlinear system R{x, u} of PDEs (II. ip with highest derivative 
of order k, solve the determining equations (13.21) for admitted sets of multipliers {A„[{7]} 
depending on the independent variables x = (x\, . . . ,x n ), the dependent variables U = 
(U 1 , . . . , U m ) and their first derivatives d x U = (U^., . . . , U™) . Two cases arise, depending 
on whether R{x, u} admits any set of multipliers of the form (13. lOj) with functions v = 
(f 1 (X), . . . ,v M (X)) satisfying some linear system of PDEs of the form (13. lip in terms of 
specific independent variables X = (Xi(x, U), . . . , X n (x, U)). 

Case I. No set of multipliers of the required form are admitted. Then one concludes from 
Theorem 3.1 that R{x, u} cannot be mapped invertibly by any point transformation to a 
linear system. 

Case II. There is an infinite set of multipliers with the required form. Typically in this 
case, the variables X,(x, U) are found directly from the specific form of the multiplier de- 
termining equations, either by inspection (see examples 4.1.2 and 4.2.2) or, more generally, 
through integration of characteristic first-order linear PDEs contained in the system of mul- 
tiplier determining equations (see example 4.3.2). Then from the admitted multipliers (13.101) 
one reads off (by Theorem 3.1) 

(a) the independent variable part of the linearizing point transformation (I3.15p . i.e. Xi —>■ 

zt = Xi(x,u); 

(b) the specific linear target system which is given by the adjoint of the linear system 

(GDH). 

Step 2: Assuming that the necessary conditions of Theorem 3.1 hold (Case II), one 
proceeds by replacing v (X) by an arbitrary function V(X) in the conservation law (13.11) 
arising from the multipliers (I3.10p for the nonlinear system H{x,u} of PDEs (II. ip . This 
will automatically yield — after integration by parts if necessary (with the use of the trans- 
formation formula for divergence expressions derived in [16]) — a conservation law identity 
of the form (13.141) for the augmented system consisting of the given nonlinear system (11.11) 
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and the linear system (13. lip . From the components of the multipliers for the linear system 
(13. lip in this identity, one directly obtains 

(c) the dependent variable part u a — > w a = W a (x, u) of the linearization mapping to the 
linear target system. 

In summary, when there exists a linearization of a given nonlinear PDE system R{x, u] 
(m > 1) by an invertible point transformation then the necessary conditions stated in Theo- 
rem 3.1 yield a linear target system along with the independent variables of this system, while 
the sufficient conditions stated in Theorem 3.2 yield the dependent variables of the linear tar- 
get system, which completes the transformation. (Examples are illustrated in sections 4.1.2 
and 4.3.2.) Furthermore, if the necessary conditions do not hold then no linearizing point 
transformation exists. 

It is easy to extend these results to include the linearization of scalar PDEs (m — 1) 
by contact transformations, where X, and W a can now depend on first derivatives of U. 
Since in this case it follows that the Jacobian J[U] and the factors may depend on 

second derivatives of U, one must consequently seek sets of multipliers {A„[£/]} depending 
on x, U, d x U and d\U . Then in Step 1 any admitted multipliers of the form (I3.10p will 
yield the independent variable part of a linearizing contact transformation, 
Xi(x, u, d x u); similarly, in Step 2 the resulting conservation law identity (13.141) will yield the 
dependent variable part of this mapping, i.e., u a —>■ w a = W a (x,u,d x u). (An example is 
illustrated in section 4.2.2.) If no such multipliers are admitted then no linearizing contact 
transformation exists. 

We remark that in certain special cases a lower differential order is possible in the required 
form of conservation law multipliers (I3.10p for the existence of a linearizing transformation. 
Such situations are important computationally because the multiplier determining equations 
will then become more over-determined and hence easier to solve. 

(1) For point or contact transformations that involve no change in the independent vari- 
ables, = Xi, one can easily show from the results in [16] that the Jacobian 
lETS]! will be trivial, i.e. J = det(M;/ Dxj) = det(5{) = 1, while the factors Q£[U] in (jOl 
will have no dependence on d x U in the multicomponent case (m > 2) and no dependence on 
d 2 JJ in the scalar case (m = 1). Hence, in this situation, conservation law multipliers can be 
sought in the form 

K[U] = 5 x y(x)QZ[U] 

depending at most on x, U and d x U if m = 1 or only on x and U if m > 2. (A fully 
algorithmic version of the linearization steps applied to this situation was given in [10], 
which allows one to detect and explicitly find all linearizing mappings that involve only a 
transformation of dependent variables.) See example 4.1.2. 

(2) Suppose a given nonlinear PDE system (II. ip of order k admits integrating factors 
{\%(x,U)} for its highest derivative terms, so that the system has the form A^ff/JG^ff/] = 
DGt[U]/ Dxi + Go[U] where Gq and G\ {% — 1, . . . , n) are of order at most A; — 1 in derivatives 
of U. Then it is straightforward to show that the factors relating the nonlinear and linear 
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PDE systems in (13.71) will be given by Q„[U] = \%(x, U) J[U] 1 , and consequently one can 
seek conservation law multipliers of the form 

K[U] =6 Xl tV X (X)\Z(x,U) 

where -Xt[£/] depends at most on x, U and d x U if m — 1 or only on x and U if m > 2 
(corresponding to linearizing contact transformations or point transformations, respectively). 
See examples 4.2.2 and 4.3.2. 



4 Examples 

To illustrate the conservation law approach for obtaining linearizations and contrast it with 
the symmetry approach, we now consider three examples. 

4.1 Linearization of Burgers' Equation 

Our first example is a nonlinear PDE system with independent variables (x\,X2) = (x,t) 
and dependent variables (m 1 ,^ 2 ) as given by 

G>] = ^_ 2Ml = , G>] = ^_ 2 ^ + ( „ lf = , (4 ,) 

where u 1 satisfies Burgers' equation 

ul x -u l ul-u\=Q. (4.2) 

4.1.1 Linearization Through Admitted Point Symmetries 

By a direct computation (solving the symmetry determining equations), one finds that equa- 
tion (14. 2p admits at most a finite number of contact symmetries. Hence there does not exist 
any contact (or point) transformation that linearizes Burgers' equation. 

On the other hand the nonlinear system (14.11) is found to admit an infinite set of point 
symmetries given by the infinitesimal generator [17,18] 

X = e u2 ' 4 (W(x, t) + g(x, ty)^ + Ag(x, t)-^j (4.3) 

where g(x,t) is an arbitrary solution of the linear PDE 

9xx = 9t- (4.4) 

Consequently, one can apply Theorems 2.1 and 2.2 to linearize the system (14.11) as follows. 
First we see from a comparison of (14. 3 p with expressions (12. 6ft that one has F 1 = g x (x,t), 
F 2 = g(x,t), with a {j = 0, (3\ = 2e" 2 / 4 , (5\ = u l e u2 ' A } f3 2 = 0, /? 2 2 = Ae u2 ' A ] the associated 
linear homogeneous system (12. 8p has X\ = x, X 2 = t as functionally independent solutions 
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and the corresponding linear inhomogeneous system (I2.9P has a particular solution (ip 1 ,^ 

/I 1 _*, 2 ! a _n,^lA\ n-n r A— I TT i , t , t~i / n1 nO\ , • r> ,1 T 



2^ 



■2 



-e u / 4 ). Then from (14.41) we have that F = {F l , F 2 ) satisfies the linear system 



F 2 X =F\ Fl = Ff. (4.5) 
Hence one obtains the invertible point transformation 

z x = x, z 2 = t, w 1 = i^e"" 2 / 4 , w 2 = -e-" 2 / 4 (4.6) 

mapping the given nonlinear system (14.11) into the linear system 

dw 2 ! dw 1 dw 2 
dx dx dt 

where both w 1 and w 2 consequently satisfy the heat equation 
d 2 w 1 dw 1 d 2 w 2 dw 2 



dx 2 dt °' dx 2 dt °' (4 ' 8) 
We thus note that the mapping (14. 6 j) and the linear system (14.71) together yield 

i 2 dw 2 /A , 

w 2 ox 

which is the well-known Hopf-Cole transformation that maps any solution of the heat equa- 
tion (14 .8p into a solution of Burgers' equation (14. 2\\ . 

4.1.2 Linearization Through Admitted Conservation Law Multipliers 

A straightforward computation of conservation law multipliers {Ai(x , t , U) , A 2 (x , t , U)} 
shows that the nonlinear system (14.11) admits an infinite set of multipliers given by 

Al = \U l e- u2 "f{x, t) + e~ u2 ^f x (x, t), A 2 = e~ u2 ^f(x, t), (4.10) 

where f(x,t) is an arbitrary solution of the linear scalar PDE f xx + f t = (namely, the 
backward heat equation). Hence, the necessary conditions of Theorem 3.1 for the existence 
of an invertible mapping of the nonlinear system (14.11) to a linear target system are satisfied, 
where, through a comparison of (I4.10p with (I3.10p . the adjoint target system is given by 

dv 1 v 2 _ q d v<2 _|_ ® vl _ q ^4 n x 

dx ' dx dt 

with dependent variables 

v = (v\v 2 ) = (/,/,) (4.12) 
and the same independent variables (x,t) as the given system (14. ip . 
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In the conservation law (13. ip arising from the multipliers (14.101) . (14.111) and (14.121) . we re- 
place {y 1 (x, t), v 2 (x, £)) by arbitrary functions (V Al (x, t), V 2 (x, £)). This leads to the following 
conservation law identity for the augmented system consisting of the given nonlinear system 
(14. ip and the linear system (14. lip : 



(y\\U l e- lj2 ^) + V 2 e-^/ 4 ) G^U] + V l e~ u '' </4 G 2 [U] 



-^v.(£-v)-^(£ + £) 

= #- (e- u2 /\-W 2 - 2U 1 V 1 )) + # (UyV^ 4 ) . 
Dx V J Dt \ J 

Consequently, directly comparing (14.131) with the identity (I3.14p in Theorem 3.2, we obtain 
an invertible mapping (m 1 ,^ 2 ) — > (w 1 = 2u 1 e~ u ^ 4 ,w 2 = 4e~" / 4 ) of the nonlinear system 
(14. ip to a linear target system given by the adjoint of the linear system (14. lip . In particular, 
this yields the point transformation (14.61) that maps the nonlinear system (14.11) into the linear 
system (14.71) up to a suitable rescaling of the dependent variables w = {w 1 ,^) 2 ). 

4.2 Linearization of a pipeline flow equation 

As a second example, we consider the pipeline flow equation [12] 

G[u] = u t u xx + u x p = 0, p = const. (4-14) 

4.2.1 Linearization Through Admitted Contact Symmetries 



The nonlinear scalar PDE (I4.14p admits an infinite set of contact symmetries whose infinites- 
imal generator is given by 

X=— ——+ (f-u — ^ — + — — (4 15) 

du x dx \ x du x ) du dt du t 

where F = F(t,u x ) is an arbitrary solution of the second-order linear PDE 
d 2 F dF 

^-m-°- < 4 - 16 > 

Consequently, one can linearize the PDE (14.141) by applying the generalization of Theo- 
rems 2.1 and 2.2 for seeking a linearizing contact transformation. In particular, X\ = u x , 
X 2 = t are found to be functionally independent solutions of the appropriate generalization 
(cf. [1,2]) of the linear homogeneous system (12.81) . while the corresponding generalization 
of the linear inhomogeneous system (12.91) has the particular solution ip = u — xu x (from 
which one deduces (p x ,p 2 ) = (u t , —x) through the contact condition in Theorem 1.2). Con- 
sequently, one obtains an invertible mapping /i given by the contact transformation 

zi = t, z 2 =u x , w = u- xu x , w Zl = u t , w Z2 = -x, (4.17) 
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which maps the nonlinear PDE (14.141) to the linear PDE given by 



4.2.2 Linearization Through Admitted Conservation Law Multipliers 

Through a computation of conservation law multipliers A(x,t, U, U x , U t ), one finds that the 
given nonlinear PDE (14.141) admits an infinite set of multipliers 

A = v(U x ,t), (4.19) 

where this function satisfies the linear scalar PDE 

v t + (U x p v) UxUx = 0. 

By inspection, since v depends on two variables (i.e., the same number of independent 
variables as U has), we see that the necessary conditions of Theorem 3.1 for the existence 
of an invertible contact transformation that linearizes the nonlinear scalar PDE (14.141) are 
satisfied with the adjoint target system given by 

^ + ^1 = (4 20) 

dr + dx 2 u 1 Uj 

for v = v(X,T) in terms of the independent variables 

X = U X , T = t. (4.21) 

In the conservation law (13.11) arising from the multipliers (14.191) . we replace v(X,T) by 
an arbitrary function V(X,T). This leads to the following conservation law identity for the 
augmented system consisting of the given nonlinear PDE (I4.14p and the linear PDE (14.201) : 

fdV d 2 (X p V) 
VG[U]-(xU x -U)J[U] — + 



^dT OX 2 

■§^((xU x - U)(U tx V + UIV X ) + ({l-p)xU x + P U)ur l V) + ^- t (lf xx (U - xU x )V 

(4.22) 



where 

D(X,T) 



J[U\ 



D(x,t) 



U x 



is the Jacobian of ( 14.211) . (For verifying the identity ( I4.22p . we note that V x = VxU xx , 
V t = VxU xt + Vt-) Consequently, the sufficiency conditions of Theorem 3.2 hold for the 
existence of an invertible mapping by a contact transformation of the nonlinear PDE (14.141) 
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to a linear target PDE which is the adjoint of the linear PDE (I4.20p . In particular, from 
a comparison of (I4.22p with the identity (I3.14p . it follows that the contact transformation 
determined by 

zi = X = u x , z 2 = T = t, w = xu x — u (4.23) 
maps the nonlinear pipeline equation (I4.14p invertibly into the linear PDE 

d 2 w dw 

x 'a>c- -W = °- < 4 ' 24 > 

4.3 Linearization of a Nonlinear Telegraph Equation 

For a final example, we consider a nonlinear telegraph system given by 

G >l = ¥-t^°' *M = £ + «V -1) = (4,5) 

with dependent variables (u l ,u 2 ) and independent variables (xi,x 2 ) = {x,t). Note that u 1 
satisfies the nonlinear telegraph equation 

(u 1 )- 2 — + l-(u 1 )- 1 =0. 



dx 2 dt \ dt 
4.3.1 Linearization Through Admitted Point Symmetries 

The nonlinear system (14.251) admits an infinite set of point symmetries given by the infinites- 
imal generator 

X = F\X, T)^- + e-tF 2 (X, T)| + e^F 2 {X, T) A + F\X, T) ± (4.26) 

where 

X I = X = x-u 2 , X 2 = T = t-logu\ (4.27) 
and (F 1 (X, T), F 2 (X, T)) is an arbitrary solution of the linear system of PDEs 

ar dx ' ox sr 1 ' 

Consequently, one can apply Theorems 2.1 and 2.2 to linearize (I4.25p . In particular, we 
see that by comparing (I4.26P with expressions (12. 6p . one has an = /3 2 = 1, au = «2i = 
Pi = Pi = 0; a 22 = e _i ,/?2 = e~ l u] the associated linear homogeneous system (12. 7p has 
Xi = x — u 2 ,X 2 = t — logw 1 as functionally independent solutions and the corresponding 
linear inhomogeneous system (I2.8P has a particular solution (ip 1 ,^ 2 ) = (x,e l ). Hence the 
mapping given by 

z\ = x — u 2 , z 2 =t — \ogu 1 , w 1 = x, w 2 = e* (4.29) 
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is an invertible point transformation under which the nonlinear telegraph system (14.251) is 
mapped to the linear PDE system given by 

^-e*^ = 0, ^!-e-^ = 0. (4.30) 

OZ2 OZ\ OZ\ oz 2 

4.3.2 Linearization Through Admitted Conservation Law Multipliers 

Through a computation of conservation law multipliers of the form {Ai(x, t, U 1 , U 2 ), 
A 2 (x, t, U 1 , U 2 )}, one can show that the nonlinear telegraph system of PDEs (I4.25P ad- 
mits an infinite set of multipliers yielding a linearization as follows. After some integrability 
analysis of the multiplier determining equations (e.g. using Crack [11]), one first obtains 

Ai = A/ 2 , A 2 = fm, (4-31) 
in terms of a function f(x, t, U 1 , U 2 ) satisfying the three linear PDEs 

fx + fu* = 0, f t + U l f m = 0, (4.32) 

(U^fmm + 2C/ X /^ - W = 0. (4.33) 

To proceed we integrate the pair of first-order PDEs ( 14.321) . which yields a reduction of the 
independent variables such that / = f(X,T) where 

X = x-U 2 , T = t-\ogU 1 . (4.34) 

Then the second-order PDE (I4.33P combined with (I4.3ip gives the infinite set of multipliers 

Ai = -fx{X,T), A 2 = -f T (X,T)/U\ fxx-f T T + f T = 0. (4.35) 

Finally, comparing (14.351) with (I3.10P in Theorem 3.1, we see that the necessary conditions 
for the existence of an invertible mapping of the nonlinear system (14.251) to a linear target 
system are satisfied, with the adjoint target system being given by 

v = (v\v 2 ) = (-f x ,-f T ). (4.36) 

and 

dv 1 dv 2 9 dv 2 dv 1 

W~W + V ax-sr=°- < 4 - 37 > 

In the conservation law (13.11) arising from the multipliers (14.351) . (I4.36P and (I4.37p . we 
replace (v 1 (X,T),v 2 (X,T)) by arbitrary functions (V l (X, T), V 2 (X, T)). This leads to the 
following conservation law identity for the augmented system consisting of the given nonlinear 
system (14. 25ft and the linear system ( 14. 3 7ft : 

v^[u] + v 2 (ur^[u} - u^j[u] + *m (j^ - ^) 

D ( rrlf 9U2 9Ul rrl\ „2 ( , T rU 1 ^ rrl 9 ^' 
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(4.38) 



where, from (14.34j) . 

D(X,T) 



J[U] 



D(x,t) 



(U 



In— 1 



dU 2 
dx 



U 1 



dU 1 
~dt 



dU 2 dU 1 
dt dx 



Consequently, by comparing f!4.38j) with the identity (13.141) in Theorem 3.2, we obtain an 
invertible mapping of the nonlinear system (14.251) to a linear target system which is given 
by the adjoint of the linear system (I4.37j) . In particular, the point transformation 



u 



Z\ = X = x — u 2 , z 2 = T = t — logw 1 , w 1 = x, w 2 
maps the nonlinear telegraph system (14.251) invertibly into the linear system 



dw 1 dw 2 



dX dT 



— w 



0. 



dw 2 dw 1 



dX dT 



0. 



(4.39) 



(4.40) 



Note that the further point transformation w 1 = w 1 , w 2 = e T w 2 maps this linear system 
(I4.40p into the equivalent linear system (I4.30p obtained from linearization of the nonlinear 
telegraph system (I4.25P through its admitted infinite set of point symmetries (I4.26p . 



5 Concluding Remarks 

An important question is which of the two methods presented in this paper is better for 
linearization? It would appear likely that the approach using admitted conservation law 
multipliers will be better computationally since, in general, the solution space for the deter- 
mining equations for conservation law multipliers is smaller than that for point (or contact) 
symmetries. A way of seeing this is to consider the case when a given system of PDEs is vari- 
ational, i.e., its linearization operator (Frechet derivative) is self-adjoint, so then Noether's 
theorem is applicable. For any such system its conservation law multipliers correspond to 
a variational subset of its admitted point (or contact) symmetries and hence, in general, 
the determining equations for conservation law multipliers are then more over-determined 
than those for point (or contact) symmetries since they include the symmetry determining 
equations subset [5-7]. 

An interesting computational way for checking for the possibility of linearization of a 
given system of nonlinear PDEs, especially for classification problems, is to apply the work 
in [19,20] to find the size of the solution space of a given system of determining equations 
without the need for obtaining any of its actual solutions. In particular, if the size of 
the solution space of the determining systems either for conservation law multipliers or for 
point (or contact) symmetries is finite dimensional, then no linearization by an invertible 
mapping is possible. But the converse (that the solution space is infinite dimensional) is 
not sufficient to imply the existence of a linearization by a point or contact transformation. 
Specifically, the solution space must have a sufficiently large number of parameters such that 
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the number of functions and independent variables, as well as the number of linear PDEs, 
matches the corresponding number in the given nonlinear system of PDEs. This counting can 
be performed fully algorithmically [19,20] provided a differential Grobner basis is available. 
However, it is important to recognize that in practice the computational complexity of finding 
a differential Grobner basis will be very high. 

It would also be interesting to develop a fully algorithmic version of the two linearization 
methods themselves. This would entail, firstly, the use of a differential Grobner basis compu- 
tation to reduce the determining equations into a consistent over-determined linear system 
satisfying all integrability conditions and having a sufficiently large solution space. Secondly, 
a strictly computational procedure must be specified to extract from this linear system the 
form of the linearizing transformation for the dependent and independent variables in the 
given nonlinear system of PDEs. 

The algorithmic steps we have used in our present work (e.g. for linearizing Burgers' 
equation, the pipeline equation, and the nonlinear telegraph system) are less automatic than 
finding a differential Grobner basis but have much less computational complexity in practice. 
In particular, they take advantage of the integration capabilities of Crack [11], including 
algorithms for direct and indirect separation of equations, for dropping redundant functions, 
for checking integrability conditions, and others. Furthermore, Crack also handles solving 
the first-order linear PDEs as needed in finding the transformations of the independent 
variables in the linearization mapping. 

Finally, we mention that a hodograph-type approach for linearization has been given 
some time ago in work [21] of Varley and Seymour on linearizing the nonlinear telegraph 
system (example 4.3). Their clever procedure for finding a linearization is not symmetry- 
based; it only applies to specific types of nonlinear PDEs and cannot determine, in general, 
whether a given system of nonlinear PDEs admits a linearization by an invertible mapping. 
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